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Abstract 

Motivated by recent developments in the AdS/CFT correspondence, we provide sev- 
eral alternative bulk descriptions of an arbitrary Wilson loop operator in Chern-Simons 
theory. Wilson loop operators in Chern-Simons theory can be given a description in 
terms of a configuration of branes or alternatively anti-branes in the resolved conifold 
geometry. The representation of the Wilson loop is encoded in the holonomy of the gauge 
field living on the dual brane configuration. By letting the branes undergo a new type of 
geometric transition, we argue that each Wilson loop operator can also be described by a 
bubbling Calabi-Yau geometry, whose topology encodes the representation of the Wilson 
loop. These Calabi-Yau manifolds provide a novel representation of knot invariants. For 
the unknot we confirm these identifications to all orders in the genus expansion. 
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1 Introduction, summary and conclusion 

The study of quantum gravity using a holographic formulation requires understanding 
how the bulk gravitational physics is described by the boundary gauge theory. Given the 
ubiquitous nature of Wilson loop operators in gauge theory, it is of interest to understand 
how Wilson loop operators in a gauge theory participating in the holographic formulation 
of a bulk theory are encoded in the corresponding bulk description. Since gauge theories 
can be formulated using Wilson loops as the basic variables, an improved understanding 
of holography may arise from identifying these variables in the bulk theory. 

This program has recently been completed in pQ for all half-BPS Wilson loop operators 
in M = 4 SYM (for previous work see [2J, [3] ) . This extends the bulk description of a 
Wilson loop in the fundamental representation in terms of a string world-sheet [U [5] to all 
other representations. It was found that for each Wilson loop, there is a bulk description 
either in terms of a configuration of D5-branes or alternatively in terms of a configuration 
of D3-branes in AdSsxS 5 (see [T] for details). The equations determining the supergravity 
background produced by these D-branes were found in [3l E] (see also [7j for the closely 
related equations for half-BPS domain wall operators), generalizing the supergravity 
solutions of LLM [8] dual to half-BPS local operators to half-BPS Wilson loop operators. 
All these asymptotically AdSsX S 5 solutions encode in their nontrivial topology the 
information about the dual operator. For other work see e.g. [9], [TOj EH [121 131 El [15] - 

A very interesting and calculable example of a holographic correspondence is the one 
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found by Gopakumar and Vafa jj6] in topological string theory 3 ]. The duality states that 
the A- model topological string theory on the deformed conifold T*S 3 in the presence of 
N D-branes wrapping S 3 is dual to the A-model topological string theory on the resolved 
conifold geometry, with the complexified Kahler modulus given by t = g s N, where g s is 
the string coupling constant on both sides of the correspondence. As shown by Witten 
[T9] , the physics on the deformed conifold is described by U (N) Chern-Simons theory on 
S 3 . 

The goal of this paper is to give the bulk description of Wilson loop operators in 
Chern-Simons theory^. As we shall see, a picture closely related to the one for M = 4 
SYM emerges, with the added benefit that we have control over all quantum corrections 
for the case of Chern-Simons theory. We find that Wilson loops in Chern-Simons can be 
described either in terms of a configuration of D-branes or a configuration of anti-branes 
in the resolved conifold geometry. Moreover, by letting the branes undergo a geometric 
transition, we find that Wilson loops can be described in terms of bubbling Calabi-Yau 
geometries, with rich topology. 

We first identify the holographic description of Wilson loop operators in Chern-Simons 
theory on S 3 - labeled by a knot in S 3 and a representation of U(N) - in terms of branes 
in the resolved conifold geometry. Just like in p], we find that there are two different 
brane configurations corresponding to each Wilson loop operator. 

The information about the knot a C S 3 is encoded in the choice of a Lagrangian 
submanifold L in the resolved conifold geometry, as explained by Ooguri and Vafa in [20] 
(see also e.g. [2TJ [22} [23]). L ends on the knot a on the S 3 at asymptotic infinity of the 
resolved conifold geometry, where the holographic dual Chern-Simons theory lives^. 

We show that the information about the representation R of the Wilson loop operator 
- given by a Young tableau (see Figured]) - is encoded in the eigenvalues^ of the holonomy 
matrix obtained by integrating the gauge connection on the branes wrapping L - denoted 
by A - around the non-contractible /?-cycle3 of the Lagrangian submanifold L. 

We show that a Wilson loop labeled by a knot a and a representation R given by 
Figure [His described either by a configuration of M D-branes (L Xl , L X2 , . . . , L XM ) or by a 
configuration of P anti-branes (L yi , L y2 , . . . , L yp ), where L Xi /L y . denotes a D-brane/anti- 

3 See [17| for an overview of results and the book [18] by Marino for a comprehensive introduction. 
4 As suggested in [IB], the Wilson loops can be described by string world-sheets in the resolved conifold. 
We elaborate on this description in Appendix [XI 

5 Taubes |24j proposed a procedure to construct the Lagrangian submanifolds corresponding to the 
knots that are invariant under the antipodal map of S 3 . The construction of the Lagrangians corre- 
sponding to arbitrary knots was achieved by Koshkin |25j . 

6 Geometrically, these eigenvalues correspond to the positions of the branes up to Hamiltonian de- 
formations (see Appendix [CJ . 

7 L has the topology of a solid torus with a boundary at infinity given by a T 2 with canonical a and 
(3 one cycles. The (3 cycle is non-contractible in L while a may or may not be contractible. 
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Ri 
R 2 



R, 



Figure 1: A Young tableau with P < N rows and M columns labeling a representation 
of U(N). Ri is the number of boxes in the i-th row and satisfies Ri > Ri+\. R T is the 
tableau conjugate to R, obtained by exchanging rows with columns. 



brane with world-volume L and with a non-trivial holonomy, shifted by the integral of 
the Kahler form J , 



X% — 



y, 



P=dD 



A,-\ l o J = g.\Ri -i + M + -J, i = l,...,M, (1.1) 
A,-\ I J = gs[Ri-i + P + l) , i = l,...,P, (1.2) 



lp=dD JD \ 

respectiveljo 

For the case of the simplest knot - the unknot - we explicitly check this identification. 
We show that this Wilson loop operator in a representation given by Figured] corresponds 
to the configuration of M D-branes (L X1 , L X2 , . . . , L XM ) 



xm x 2 X\ 



Figure 2: Brane configuration in resolved conifold describing Wilson loop in a represen- 
tation given by Figure CD 

where L Xi denotes a D-brane wrapping a Lagrangian submanifold L sitting at the position 
Xi = g s (Rf — i + M + 1/2) on the inner edge of the toric diagram of the resolved conifold. 
In particular, a single D-brane on the inner edge corresponds to a Wilson loop in the 
antisymmetric representation. Since the effective size of the inner edge of the resolved 

8 l 7 is integrated over a disk D. D is a relative cycle in the resolved conifold and has (3 C L as its 
boundary. 
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conifold is t e ff = g s (N + 1), we recover (from x\ < t e ff) the group theory bound R\ < N 
for the number of boxes in a column from the compactness of the P 1 of the resolved 
conifold. 

The same Wilson loop operator can also be described by a configuration of P anti- 



Figure 3: Brane configuration in resolved conifold describing Wilson loop in a represen- 
tation given by Figure 1. 

where L y . denotes an anti-brane wrapping a Lagrangian submanifold L sitting at the 
position yi = g s (Ri — i + P+1/2) on an outer edge of the toric diagram. In particular, a 
single D-brane on the outer edge corresponds to a Wilson loop in the symmetric repre- 
sentation. Since the outer edge is non-compact, R\ can be an arbitrarily large integer, as 
expected from group theory. The bound on the number of anti-branes P < N can also 
be understood geometrically. Once we put P anti-branes in the resolved conifold, the 
effective size of the space is t e ff = g s (N — P), thus recovering the group theory bound 



These results are obtained by computing the A-model topological string partition func- 
tion in the presence of these branes and showing that it agrees with the Chern-Simon 
results of Witten [19]. In the crystal melting description [261 [27] of these amplitudes, the 
Young tableau of the dual Wilson loop has a nice geometrical realization in the crystal. 
The corners produced in the crystal by the insertion of branes [28l [27] has exactly the 
shape of the corresponding Young tableau! 

We also find a purely geometric description of the unknot Wilson loop operators in 
terms of Calabi-Yau geometries that are asymptotic to the resolved conifold. We show 
that the backreaction of the brane configurations in Figures 2 and 3 can be taken into 
account exactly, giving rise to a Calabi-Yau geometry without D-branes! The two dif- 
ferent D-brane configurations in Figures 2 and 3 yield the same Calabi-Yau geometry, 
and the information about the representation R of the Wilson loop is now encoded in 
the topology of the corresponding Calabi-Yau. Physically, the appearance of this rich 
class of Calabi-Yau geometries can be understood as arising from a new class of geo- 
metric transitions, where non-compact Lagrangian D-branes/anti-branes are replaced by 
non-trivial topologies supporting the "flux" produced by the branes. The computation 



branes (L yi , L 





P< N. 
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of the A-model topological string amplitude on these "bubbling" Calabi-Yau geometries, 
obtained by letting the branes undergo a geometric transition, exactly reproduces the 
result for the expectation value of the corresponding Wilson loop obtained by Witten 
[T9] . Therefore, these bubbling Calabi-Yau geometries provide a novel representation of 
knot invariants of Chern-Simons theory. 

The topology of the Calabi-Yau geometry corresponding to a Wilson loop is best 
understood by giving the following parametrization of a Young tableau^ 




Figure 4: The correspondence of the Wilson loop in representation R and the bubbling 
geometry. In (a) the Young tableau R, shown rotated, is specified by the lengths k of 
all the edges. lim+\ is N minus the number of rows. Equivalently, l{ are the lengths 
of black and white regions in the Maya diagram. The Wilson loop in representation R 
(a) is equivalent to the toric Calabi-Yau manifold given by the web diagram (b). There 
are a total of 2m + 1 bubbles of P 1 in the geometry. The sizes of the P^s are given by 
t { = g s k, i = 1, ...,2m + 1. 

The Calabi-Yau geometry dual to a Wilson loop labeled by a Young tableau given in 
Figure H^a) is given by the toric web diagram shown in Figure H^b). The corresponding 
Calabi-Yau geometry has 2m + 1 nontrivial P^s. The associated complexified Kahler 
moduli are given by ti = g.Ji for i = 1, . . . , 2m + 1, where U are the integers in Figure 
IHa) parametrizing the Young tableau. Similarity between bubbling AdS geometries and 
toric geometries was noted in [8]. 

Although the explicit quantitative checks we make are all based on the unknot, our 
physical arguments apply to arbitrary knots and links. Indeed it is possible to write 
the closed Gopakumar-Vafa (GV) invariants [IE] associated to the bubbling Calabi-Yau 
corresponding to an arbitrary knot in terms of the open GV invariants [20| [22] associated 

9 Informally, l oc id is the number of rows in the tableau with the same number of boxes while l even is 
the number of columns in the tableau with the same number of boxes. 
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to the Lagrangian branes giving rise to the bubbling Calabi-Yau after geometric transi- 
tion. We plan to report on this relation and on more general geometric transitions in the 
forthcoming paper [29J . 




The extent of the analogy between AdS / CFT and the GV duality is striking. We hope 
that the very explicit realization of the holographic correspondence found in this paper 
- including all quantum corrections - provides us with hints for deeper understanding of 
the AdS/CFT correspondence and holography in further detail. 

The plan of the rest of the paper is as follows. In section 2 we integrate out the degrees 
of freedom introduced by a configuration of non-compact branes intersecting S 3 along a 
knot a and show that the net effect is to insert a Wilson loop operator for U(N) Chern- 
Simons theory on S 3 . In section 3 we identify the brane configurations in the resolved 
conifold dual to a Wilson loop operator. The identification is explicitly checked for the 
case of the unknot. In section 4 we show that the Wilson loop operators can be related 
to the closed string partition function of bubbling Calabi-Yau's, which are interpreted 
as arising via a geometric transition of the brane configuration in section 3. Appendix 
lAl revisits the description of Wilson loops in terms of string world-sheets. Appendix 
iBl reviews the melting crystal models found in [27] . The explicit computations for the 
unknot, which support our proposal in this paper, are performed using the crystal mod- 
els. Appendix O explains relevant details of the target space theory of open topological 
strings. 

2 Wilson loops as branes in the deformed conifold 

In this section we find a brane configuration in T*S 3 whose effective description is given by 
a Wilson loop operator of Chern-Simons theory on S 3 . The basic idea is to add extra non- 
compact branes beyond the N D-branes wrapping the S 3 that support the U(N) Chern- 
Simons theory participating in the large N duality. We show that integrating out the 
degrees of freedom introduced by the extra branes has the effect of inserting into the U(N) 
Chern-Simons path integral on S 3 a Wilson loop operator in a particular representation 
R of U (N) determined by data of the brane configuration. Having identified the brane 
configuration in T*S 3 corresponding to a Wilson loop, in the next section we find the 
holographic description of Wilson loops in terms of branes in the resolved conifold. 

A Wilson loop operator of Chern-Simon theory on S 3 is labeled by a representation R 
of U(N) and by an embedded oriented circle a - also known as a knot - in S 3 . It is given 




by 




(2.3) 



where P is the path ordering operator. 
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A microscopic description of a Wilson loop operator can be given by adding new degrees 
of freedom localized on the knot on which the operator is supported. In the string theory 
realization, the new localized degrees of freedom arise microscopically by quantizing open 
strings between the D-branes wrapping S 3 and new branes intersecting the S 3 along the 
knot a. 

As shown by [20], a Lagrangian submanifold L in T*S 3 can be found such that it 
intersects S 3 along an arbitrary knot a. Given a knot a in S 3 parametrized by q l {s), 
where q l are coordinates on S 3 , then L is determined by the following non-compact 
Lagrangian submanifolc0 

L = Pi) |^;^ = 0}, (2.4) 

i 

where Pi are coordinates in the fiber of T*S 3 . The topology of L is R 2 x S 1 or more precisely 
that of a solid torus. The brane is non-compact and has a boundary at asymptotic infinity. 
The boundary is a T 2 with canonical a and (3 one cycles. In L, the a-cycle of the T 2 is 
non-contractible while the /5-cycle of the T 2 is contractible^]. 

Therefore, we want to consider the gauge theory living on the brane configuration in 
T*S 3 given by N D-branes wrapping S 3 intersecting M branes wrapping L. The action 
is given by [20J 



where 



S = Scs(A) + Scs(A) + l X \d + A + A) X , 

J a 

Scs(A) = — f Tr ( A A dA + -A A A A A ] , 
9sJx\ 3 J 



(2.5) 



(2.6) 



and A, A are the connections on X = S 3 and X = L respectively. \ are bifundamental 
fields localized along the knot a arising from the quantization of open strings with one 
end on S 3 and the other one on L. Depending on whether we wrap a D-brane or an 
anti-brana^ 2 ! on L, \ is either a fermionic or bosonic field. 



10 We recall that the Kahler form of T*S 3 is given by uj = £? =1 dp t A dq\ 

11 As we will see in the next section, the roles of a and (3 are reversed in the corresponding Lagrangian 
submanifold L in the resolved conifold; (3 being non-contractible in the corresponding L. 

12 In topological string theory an anti-brane is exactly the same [301 IS] as a ghost brane, obtained by 
reversing the sign of the D-brane boundary state. That's why the statistics of the open string fields is 
reversed when open strings are stretched between a D-brane and an anti-brane as compared to when 
they are stretched between two D-branes. 
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2.1 Wilson loops as D-branes 



If D-branes wrap L then x are fermions. Integrating out x 1S straightforward, it inserts 
the so called Ooguri-Vafa operator |20p*l 

J^TrgrPexp j> A ■ Tr Q Pexp I A, (2.7) 

Q 

where Q T is the Young tableau conjugate to Q. We now want to integrate A ou 



e iScs{A) ^2Tr QT Pexp I A- I [DA]e iScs{A) TT Q P exp I A 

r\ J a. J J a 



(2. 



In order to proceed we now briefly recall some well known facts about Chern-Simons 
theory. In the present case, Chern-Simons is defined on a solid torus with a boundary at 
infinity given by a T 2 . The path integral is a wave function for Chern-Simons theory 
on a T 2 and depends on the boundary condition imposed at infinity. 

The path integral in (12.81) has the insertion of TrgPexp^^, that is a Wilson loop 
along the non-contractible cycle a of T 2 . This creates a state of Chern-Simons theory on 
T 2 labeled by \Q). Therefore (J2HD yields 



e iScs(A) 



^Tr Q TPexp<[ A-(i>\Q). (2.9) 



J a 



The holonomies around the a and (3 cycles of T 2 are canonically conjugate [32 J. The 
path integral is defined by specifying the holonomy at infinity either around the a or (3 
cycle. Roughly speaking, the holonomy around a plays the role of the position operator 
while the holonomy around the (3 cycle plays the role of momentum. In the context of 
Chern-Simons theory, momentum can be identified with the highest weight vector of a 
representation R, shifted by the Weyl vector. 

Since our aim is to find the Wilson loop in a particular representation R f l 2 . 3 j) . we 
choose the boundary condition = (R T \. This means that we choose a boundary 
condition at infinity with non-trivial holonomy around the contractible cycle (3 in L. 
Since the states labeled by representations form an orthonormal basis, this picks out the 
Q = R T term in (12. 9p . In other words, imposing the boundary condition = (R T \ is 
equivalent to focusing on the Q = R T term in (12.91) . 

We now give the physical interpretation of this boundary condition. In a nutshell, the 
holonomy around the contractible cycle (3 measures fundamental string charge. 
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The path integral over x reduces to that of free fermions. This formula can then be derived by 
noting that fl/lil 1 + x iUi) = Eq TtqtXTtqY, where X = diag(a; l ), Y = diag(yj). 

14 In this analysis we omit the path integral over A, which is to be done at the end. In [20] the path 
integral over A was performed while the integral over A wasn't. 
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Open String Endpoints as Anyons 



To understand the physical meaning of the boundary condition, let us consider the 
effect of the Wilson loop Tr^ exp § A on the non-compact branes. 

Consider first the case when M — 1, i.e. for a single brane on L. The action of the 
brane in the presence of the Wilson loop (swept by an open string end point) is given by 



S = Scs(A) +kf A. (2.10) 

Therefore, the corresponding equations of motion are given by 

F z - Z = g s k5 2 {z), (2.11) 

where z is the complex coordinate parameterizing the R 2 in L. Physically, the endpoint 
of an open string behaves like an anyon when viewed from the brane world-volume. 

By integrating (12.111) we conclude that having k fundamental strings ending on L 
introduces a non-trivial holonomy around the contractible cycle (3 of the T 2 in L. The 
holonomy is given by 

<pA = g s k. (2.12) 

In the general case of M arbitrary, the treatment of the non-Abelian equations of 
motion is more subtle, as noted originally by Witten [33]. Fortunately, we can borrow 
results from [32], where the insertion of the Wilson loop in the representation R T was 
found to induce the holonomy 

Ai = g s (rJ - i + X -M + , i = 1, 2, M. (2.13) 

Therefore, the boundary condition = (R T \ is tantamount to introducing the holonomy 
given by (12.131) at infinity. We denote by {L r t, . . . , L^) the configuration of branes with 
this holonomy (12.131) . 

When we go through the transition to the resolved conifold, we will see that there is a 
shift by \g s M in the effective holonomy, coming from the backreaction of the geometry 
on the branes. In Appendix C we show that the gauge invariant quantity is given by the 
holonomy shifted by the integral of the Kahler form over a disk ending on /3, thus giving 
rise to the shift b\l 15 l \g s M in the resolved conifold. 



15 The shift of the Kahler modulus due to the insertion of branes is well known [34, 28, 27] and it is 
given by g s M. Therefore, when the Kahler modulus is integrated over a disk bounding (3 - as opposed 



to over P 1 - we get a shift by \g s M. 
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Therefore, the D-brane configuration {L r t, . . . , L r t ) we have considered in T*S 3 in- 
serts a Wilson loop operator along the knot a 

e iS C s(A) . Trfi p exp I a. (2.14) 

J a 

Summarizing, we have derived the identification 

{Ljff, . . . , Ljff) <- Tr R Pexp I A. (2.15) 



J a 



2.2 Wilson loops as anti-branes 



In topological string theory an anti-brane has the interpretation as a ghost brane [301131] . 
whose boundary state is minus that of a brane. This sign changes the statistics of 
the open string fields between branes and anti-branes. This is why the x fields in the 
brane configuration we have described in T*S 3 are bosonic when anti-branes wrap the 
Lagrangian submanifold L. 

When x are quantized as bosonic fields, integrating them out yield@ 

^Tr Q Pexp & A ■ Tr Q Pexp <t A. (2.16) 

Q J a J a 

We can now easily calculate the path integral over A following our previous discussion 
for D-branes. If we denote by (L Rl , . . . , L Rp ) a configuration of P anti-branes with 
holononr\l 17 l 



At = g s ( Ri - i + \p + i ) , i = 1, 2, .., P, (2.17) 



corresponding to the boundary condition (?/>| = (R\, we are left with the insertion of a 
Wilson loop operator along the knot a given by 

e iS cs (A) . TrflPexp I A (2.18) 

J a 

Therefore we arrive at the identification 

(L Rl , . . . ,L Rp ) ^Tr^Pexp I A. (2.19) 

J a 

We now study the bulk description of Wilson loops in Chern-Simons theory in terms 
of branes in the resolved conifold geometry. 



16 Now the path integral over \ reduces to that of free bosons. The final answer is obtained by noting 
that nLnii T^kn = Eq^qXTtqY). 

17 As before, there will be a shift by \g s P in the resolved conifold. 
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3 Wilson loops as branes in the resolved conifold 



In the previous section we have shown that a Wilson loop operator on any knot a and 
for any representation R can be obtained by integrating out the physics of a D-brane 
configuration (L r t, . . . , L r t ) or anti-brane configuration (L Rl , . . . ,Lr p ). 

To obtain the resolved conifold description of a Wilson loop we follow the brane config- 
uration through the conifold singularity. It is possible to construct a Lagrangian subman- 
ifold L explicitly for every knot a in S 3 [24l [25] . Physically, this Lagrangian submanifold 
in the resolved conifold is the dual description of the Lagrangian submanifold (12.4ft in 
T*S 3 described in the previous section. Topologically L is also R 2 x S 1 and has an 
asymptotic boundary given by a T 2 . 

The Lagrangian submanifolds constructed by Taubes have the property that they end 
on a knot a in the S 3 at infinity. From the point of view of holography, this is precisely 
as expected. The dual U(N) Chern-Simons theory lives on the S 3 at asymptotic infinity 
in the resolved conifold. Given that we are looking for the resolved conifold description 
of Wilson loops, it is expected that the bulk description is given by a bulk object which 
ends on a knot, thus introducing the appropriate source for a Wilson loop operator. 

A crucial role in the derivation in the previous section is played by the holonomies 
around the contractible cycle of the T 2 in L 

I Ai- (3.20) 

h 

This data has a nice geometrical interpretation in the resolved conifold. As one follows 
the branes through the conifold singularity, the contractible cycle (3 of the Lagrangian 
L in T*S 3 grows to become a non-contractible cycle (3 of the T 2 on the corresponding 
Lagrangian L of the resolved conifold. 

As explained in more detail in Appendix C, the holonomy of the complex gauge field 
A has the interpretation as the modulus^ of the brane. More precisely, in addition to 
the holonomy of the gauge field, the gauge invariant modulus has a contribution from 
the Kahler form J integrated over a disk ending on j3. This contributes \g s M to the 
modulus since f D J — \g s M (see footnote fT5l) . 

Therefore, the Wilson loop operator described by the brane configuration (L r t, . . . , L r t ) 
in the deformed conifold has a bulk interpretation in terms of M Lagrangian D-branes 
(L Xl , . . . ,L XM ) in the resolved conifold. Moreover, the modulus Xi of the z-th brane in 
the resolved conifold is determined by the holonomy data (I2.13P 

Xi=<f A+ [ J = g s [Rj -i + M + i j i = l,...,M (3.21) 

JdD=/3 J D V J 

18 In the toric diagram, the modulus is the position of the brane. 
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Similarly, the brane configuration (L Rl , . . . , Lr p ) in the deformed conifold has a bulk 
interpretation in terms of P Lagrangian anti-branes (L yi , . . . , L yp ) in the resolved coni- 
fold. Moreover, the modulus yi of the z-th brane in the resolved conifold is determined 
by the holonomy data (I2.17P 

Vi=(f A+ f J = g s [R l -i + P + )\ i = l,...,P (3.22) 

JdD=p Jd V 1 J 

We now proceed to verify this identification for the case when the Wilson loop operator 
is defined on the simplest knot, the unknot. When the deformed conifold geometry is 
defined by 

-21-24 - -22-23 = (3.23) 

the unknot is parameterized by (zi, z 2 , z 3 , Z4) = (0, y/JIe 1,9 , —^/Jle~ te , 0) with < 6 < 2n. 

For the case of the unknot, the Wilson loop operator Wr in (12.31) is described by the 
configuration of D-branes in Figure [jpl . The corresponding Lagrangian submanifolds 
end on the inner edge, thus accounting for the group theory bound Rj < N. This is 
consistent with the requirement that (3 should be non-contractible. 

It was shown in [27], as reviewed in Appendix [Bj that the Wilson loop vev for the 
unknot can be rewritten as 

voo e -»* /_ , . \ il ~ e ~nxi 1 e -n(t-xi) 

(W R ) = M(q)e-^^ m(l - e-^-^))l[ e xpJ2 . ( 3 - 24 ) 

\i<j J i=l n=l ^ ' 

up to unimportant factors we suppress here. Here M(q) = YYjLii^ ~ ( t')~^ is the MacMa- 
hon function, t = g s (N + M) and Xi = g s (Rf — i + M + |), i = L, M. This is exactly 
the A-model amplitude for the brane configuration in Figure 2 r°l . thus confirming our 
identification 

(W R ) = {{L X1 ,...,L XM )). (3.25) 
We also show in Appendix [B] that the Wilson loop vev can also be written as 

x -nyi _ p -n(t+yi) 



K i<j / i=l n=l ' 



19 Let us parameterize the resolved conifold by | Ai | 2 + A2 1 2 — |Ci| 2 ~ IC2I 2 = Rei up to U(l) equivalence 
with charges (1, 1,-1 — 1). In Figures [2] and [3j the left, right, top, and bottom regions correspond to 
I A x | 2 = 0, |A 2 | 2 = 0, |Ci| 2 = 0, and |C 2 | 2 = 0, respectively. L is given by |Ai| 2 -Rex = |A 2 | 2 - Re (t - x) = 
ICi| 2 — IC2I 2 , arg(£i£ 2 AiA 2 ) = n. When the conifold is singular the z and (A,C) coordinates are related 

The factor Hi<j(l ~~ e ^ x> 1S the contribution from annulus diagrams between branes. This is 
essentially the Vandermonde determinant, i.e. the Weyl denominator. It demonstrates that the branes 
are fermions, as shown first in the mirror B-model |34| . The rest of (j3.24[) , up to M (q), can be computed 
by the topological vertex technology [33] . 



as 

20 
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with yi = g s (Ri — i + P + ~), i = 1, P, t = g s (N — P). This is exactly the A-model 
amplitude for the brane configuration in Figure 3, thus confirming our identification^*] 

(W R ) = {(L yi ,...,L yp )). (3.27) 

In summary we have identified the bulk description of Wilson loop operators in Chern- 
Simons theory in terms of D-branes and anti-branes in the resolved conifold. Moreover, 
this identification has been explicitly verified for the case of the unknot. 



4 Wilson loops as bubbling Calabi-Yau's 

A geometric transition is a phenomenon in which a stack of D-branes is replaced by a 
new geometry with "flux" . More precisely, when the number of branes in the stack is 
large, the system is better described by a certain geometry where the appropriate fields 
that encode the charges of the branes are turned on. In physical string theory these fields 
are RR fluxes originally sourced by the D-branes. In topological string theory, the role 
of the fluxes is played by bulk gauge fields, namely the Kahler form and the B-field. The 
change in the geometry is such that the non-trivial cycle, which originally surrounds^! 
the branes and is homologically trivial, becomes a non-trivial cycle that supports the 
"flux". There are by now many examples of this phenomenon. Here we will find a new 
class of geometric transitions in topological string theory. 

It has been argued in 0, [6] that the D-branes realizing the half BPS straight Wilson 
lines in M = 4 SYM can undergo transition to certain bubbling geometries that asymptote 
to AdSs x S 5 . To obtain these geometries one makes an ansatz for supergravity fields 
based on the knowledge of the symmetry of the branes describing the Wilson loops. One 
then imposes the BPS condition and finds that the supergravity solution is determined by 
simple data, namely a black-and-white pattern on a line as in Figure H](a). It is expected 
that the data corresponds to the representation of the Wilson loop [3]. 

We found in section [3] that Wilson loop operators in Chern-Simons theory can be 
realized by a configuration of branes or anti-branes in the resolved conifold. When the 

21 The same amplitude would arise if the branes ended on any of the four outer edges. One can 
convince oneself that they end on the lower-left edge by the following argument. In the convention of 
footnote [Tni one can show that a is the trivial cycle in L if L ends on the lower- left or upper- right edge. 
(3 is trivial if L ends on the lower-right or upper-left edge. Since we have non-trivial holonomy of a flat 
connection along /3, [3 has to be non-trivial, and thus L has to end on either the lower-left or upper-right 
edge. By carefully keeping track of the orientation of the cycle (3, one can show that increasing the 
holonomy is gauge equivalent to moving the branes toward the lower-left. Since one can increase the 
holonomy without bound, the branes must end on the lower-left edge. Such L is given by the equations 

|Ai| 2 HA 2 | 2 -Re(t + 2/ )HCi| 2 -Rey=|C2| 2 , arg(CiC 2 A 1 A 2 ) = tt. 

22 Let W be the world-volume of the brane, and M a homologically trivial cycle. We say that M 
surrounds W if there is a chain N such that dN = M and if the intersection number of N and W is 
non-zero. 
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number of D-branes in a stack is large we expect that the system has a better description 
in terms of pure geometry. The D-branes wrap a Lagrangian submanifold L of topology 
R 2 x S 1 . The neighborhood of the submanifold, modeled by the normal bundle, is locally 
R 5 x S 1 . A contractible S 2 in the transverse R 3 surrounds the branes. The geometric 
transition of the D-branes makes the S 2 non-trivial while making the S 1 contractible. 
More precisely, the topology change is described by a surgery procedure: we cut out a 
tubular neighborhood of topology (3-ball)xR 2 x S 1 from the ambient space X and glue 
in a region of topology S 2 x R 2 x(disk) to get a new space X'. A relative 2-cycle with 
boundary on L combines with the disk to become a 2-cycle in X'. 




Figure 5: The simplest example of the correspondence between a Wilson loop and 
a bubbling geometry. The Wilson loop along the unknot in the U(N) representation 
specified by the Young tableau (a) is equivalent to the toric Calabi-Yau manifold given 
by the web diagram (b). The Kahler moduli in (b) are given by t\ = g'sJn, t 2 = gj, t 3 = 
g s (N — m). The bubbling CY manifold (b) arises from geometric transition of I branes 
in (c) as well as m anti-branes in (d). 

Let us now specialize to the Wilson loop along the unknot, in a representation of the 
form Figure 0(a). The Wilson loop is realized by a stack of I Lagrangian D-branes on 
the inner edge at positions Xj = g s (l + m — i + 1/2), i = 1, These branes live in 
the resolved conifold with Kahler modulus t = g s {N + I). We propose that when I and 
m are large, these branes undergo geometric transition to the toric Calabi-Yau manifold 
whose toric web diagram is shown in[^b). The three Kahler moduli are t\ = g s m, t 2 — 
g s l, £3 = g s (N — m). Note that the new geometry has topology that is precisely as 
described above. The contractible sphere that originally surrounded the branes become 
a non-trivial cycle of size g s l. The two holomorphic disks that originally ended on the 
branes become spheres of sizes g s m and g s (iV — m). Our proposal is supported by the 
explicit computation for the unknot. If we substitute Xi = g s (l + m — i + 1/2) into 
the open+closed string partition function (13.241) in the brane set up of Figure 0(c), it 
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become; 



oo _ 

M{q) 2 expJ2 -pj2 {~ e ~ ntl ~ e ~ nh ~ e ~ nt3 + e~ n{tl+t2) + e~ n{t2+t3) - e - n( * 1+ * 2+ ' 3) )4.28) 

n=l ^ ^ 

This is precisely the closed string amplitude for the bubbling Calabi-Yau in Figure [5(b), 
computed using mirror symmetry and integrality^! This constitutes very strong evidence 
for our proposal for the geometric transition. 

In section [3] we showed that the same Wilson loop is realized by a stack of m anti- 
branes on an outer edge with holonomies %ji = g s (l + m — i + 1/2), i — 1, ...,m turned 
on. (See Figure 13(d).) The anti-branes live in the resolved conifold with t = g s {N — m). 
These anti-branes also undergo geometric transition. The resulting geometry is the same 
as for the branes on the inner edge. In this case the 2-sphere surrounding the anti-branes 
acquires size g s m, and one 2-sphere of size gj arises from a holomorphic disk ending 
on the anti-branes. The open+closed string partition function (I3.26P for the anti-branes 
again becomes^] the closed string partition function (I4.28P after substituting the values 
for yi. 

Let us now discuss the Wilson loop along the unknot in the general representation 
R. It is convenient to parameterize R in terms of lengths U of the edges as in Figure 
|4(a). The open+closed string partition functions (I3.24p and (13.261) both become, after 
substituting the values for holonomies, 

oo / 

M(gr +1 exp^-— - £ e"^+ £ 

n=l L ' \ l<i<2m+l l<i<2m 

_ ^ ^ e -n(ti+t i + 1 +t i+2 ) _ g -n(tl + ...+*2m + l) 

l<i<2m-l 

We recognize this as the closed string partition function for the toric Calabi-Yau whose 
toric diagram is shown in Figure H(b) [31] ! This complicated geometry arises via geo- 
metric transition of branes or anti-branes. The description in terms of geometry is more 
appropriate when l iy hence the cycles, are large. These Calabi-Yau manifolds provide a 
novel representation of knot invariants in Chern-Simons theory. 

In the forthcoming paper [29] we will show that the open GV invariants for an arbitrary 
knot determine the closed GV invariants of the bubbling CY geometry as well as discuss 
generalization of the newly found class of geometric transitions. 

23 For precise agreement we should include £(<?)' in (|3.24|) where (;(q) = Yi'jLii^ ~ <7 J ) -1 - This factor 
does not contribute to perturbative amplitudes because of its modular property [28] . 
24 See subsection 9.1 of [34] 
25 Insert £,(q) m to (|3.26p for precise agreement. 
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Appendix 

A Wilson loops as world-sheet instantons 

As in AdS/CFT, Wilson loops in Chern-Simons theory are expected to be realized by 
non-compact world-sheets p21 [35] in the resolved conifold. In this appendix we explain 
how this works and point out that the description in terms of world-sheets naturally gives 
rise to the generating functional of Wilson loops originally obtained in [20J. The gener- 
ating functional is the amplitude for branes wrapping a certain Lagrangian submanifold. 
Conversely, this implies that exact computation of the brane amplitude amounts to exact 
quantization of string theory around a classical string solution. 

In Chern-Simons perturbation theory, the expectation value of the Wilson loop opera- 
tor in the fundamental representation is computed by summing Feynman diagrams with 
external legs ending on the loop: 




(A.30) 



Here the filled circles are the end points of the gauge field propagators, and are integrated 
along the Wilson loop. When the Feynman diagrams are drawn as fat graphs, they can 
directly be viewed as infinitely thin open string world-sheets [T9] . To provide for a string 
theory description, we place an extra holomorphic world-sheet that has the loop as its 
boundary because the boundary of the world-sheet is the Wilson loop. This surface is 
unique and is given as follows. The equation for the deformed conifold is given in (I3.23p . 
Let us introduce w coordinates by z% = W\ + iw2, z& = w\ — 1W2, 22 = u>3 + zu>4, z 3 = 
—W3 + 1W4. Then the equation becomes 

w x 2 + w 2 2 + w 3 2 + w A 2 = fj,, w t e C, % = 1, 2, 3, 4. (A.31) 

Let us write W{ = a* + ih{ with real Oj, hi. The S 3 is the real locus {pi = 0) of the 
hypersurface. A knot can be parameterized as 

Oi = MO), 0<9<2n, (A.32) 
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where fi(9) are real functions satisfying ^ fi(9) 2 = A*- Then the holomorphic surface is 
given by 



Thin open string world-sheets should then be glued to the boundary of the holomorphic 
surface. In other words, we are quantizing the open string theory around a classical string 
solution represented by the holomorphic world-sheet. As in AdS/CFT we need to cut off 
the infinite area of the surface. After IR regularization the contribution of the infinite 
holomorphic surface is simply e~ Area . By subtracting the area, string theory around the 
holomorphic world-sheet computes the Wilson loop vev. 

It is easy to generalize this consideration to an arbitrary multi-trace case 



k is an infinite vector of non-negative integers. If we know these expectation values we 
can compute the vev of the Wilson loop in an arbitrary representation with the help 
of the Frobenius relation. The computation of such an operator is done by expanding 
the string theory around a classical string solution in which we place k\ singly wrapped 
surfaces, k 2 doubly wrapped surfaces, etc. The classical configuration involves a total of 

kj distinct world-sheets. Note that there are quantum fluctuations that connect them. 
They represent the interactions among fundamental strings. 

The IR regularization of the world-sheet can be conveniently done by introducing M 
non-compact D-branes on which the second boundary of the large world-sheet ends. 

These D-branes wrap a Lagrangian L. It is in fact convenient to deal at once with 
configurations with all possible k (with kj = for all j > M). Let s be the area of the 
holomorphic surface. The idea is that the surface appears in the perturbative quantization 
of topological strings in the background with non-compact D-branes. Let us denote by 
V the M x M holonomy matrix for the gauge field on L along the boundary of the 
surface. Because the world- volume is non-compact, the gauge field on the non-compact 
branes is essentially non-dynamical. The contribution from the stack of world-sheet 
instantons specified by k picks up the factors rijli(TrV J ') j from the holonomy on the 
second boundary. The open string partition function Z(V) of the system with branes is 
given by 



1/zj: = ]-/ Ylj{kj\)j ki is the symmetry factor. Up to e~ s this agrees with the generating 
functional of Wilson loops obtained in [20] through consideration in the open string 
channel. Here we provided a closed string channel derivation. 



Wi(0 = fi(0, C e C, C ~ C + 2vr, ImC > 0. 



(A.33) 




(A.34) 




(A.35) 
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Let us now follow the conifold transition to the resolved side. The Gopakumar-Vafa 
duality implies that Z(V) can be identified with the partition function of topological 
strings in the resolved conifold, in the presence of non-compact branes with holonomy 
V. The contribution to the coefficient of YljLi( r ^ r ( e ~ s ^ r )'') now comes from a stack of 
world-sheet instantons. Their boundaries are on L and they are wrapped multiply in the 
way specified by k. The branes provide IR regularization. Thus our conclusion is that the 
Wilson loop YlJLi (TrPe- 7 ^" 4 ) 3 is identified with the stack of non-compact world-sheet 
instantons whose wrapping numbers are k = (ki, k^ ...). The exact computation of the 
brane amplitude (1A.35I) . which is possible for unknot, then leads to exact quantization 
of string theory around the classical string solution. Note that the interactions between 
world-sheets, corresponding to quantum fluctuations that connect world-sheets, arise in 
the perturbative computation of the brane amplitude. 



B Wilson loops as crystals 

In this section we revisit the crystal description of the Wilson loop vev for the unknot, 
developed in [27]. Besides being interesting on its own right, our aim here is to use 
such crystals to facilitate the computations showing that the Wilson loop vev equals the 
A-model amplitudes of branes, anti-branes, and the bubbling Calabi-Yau'Q. 

Let R be the representation corresponding to a Young diagram. Consider the three 
dimensional region shown in Figure [61 A unit cube in the region is considered to be 
an atom. Outside the region is the air. We regard the region as the initial crystal 
configuration, which is specified by Young diagram R in Figure Ufa). An atom can melt 
away only if three of its faces are already exposed to the air [26]. Now let us consider 
the ensemble of all crystal configurations in which some atoms have melted away. The 
partition function is defined by summing over all melted crystal configurations with 
Boltzman weight q E , where E = (number of removed boxes). It was shown in [27J that 
the unnormalized expectation value of the Wilson loop (Wr) for unknot agrees with the 
crystal partition function up to some factorial. 

This type of crystal melting problem can be treated efficiently using the free CFT 
techniques [26]. For our purposes, all the reader needs to know is the commutation 
relation of two operators T±(z) and how they act on the vacuunQ 

r + (z + )r_(z_) = - — 1 r_(z„)r + ( z+ ), r + (*)|o) = (o|r_(*) = o. (b.36) 

1 — Z + / Z- 



26 The computations can in principle be performed starting with the expression (Wr) ~ 
deti<ij<jv(g :, ' fli ~ 1 ''), though the crystals provide more transparent methods. 

27 For the prefactors that play no role for us, see [27] . 

28 If the reader wishes to understand how the crystal melting problem is translated to the free CFT 
expressions, formulas like r+(l)|i?) = J2r-<q \Q) tna t can be found in [26] . 
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Figure 6: Crystal configuration before melting. This crystal corresponds to Figure HI 



The partition function (Wp) can be expressed in several different ways. Here we con- 
sider the expressions which differ onlj@ in the way the Young diagram R is parametrized. 
The first parametrization uses the number Ri of boxes in the i-th row. The second uses 
the number Rf of boxes in the i-th column. The third uses the length li of i-the edge. 
See Figures [T] and Ufa) for examples. 

In the parametrization in terms of the number R t of boxes in the i-th row, the Wilson 
loop vev is given by 



(W R ) = (0| 



n 



xr_(g 



i=Ri 
R P -P\ 



Rp-P-X 

n r + («- 



Ri-2 

n r +(<? ? ) 

i=R 2 -l 

ffr-tf 

i=-N 



r_(g R2 - 2 )... 

|0>. 



(B.37) 



( 1B.36P allows one to compute this and obtain (13.261) that is the partition function in the 
presence of anti-branes on the outer edge. The special case P = 1 of this interpretation 
was found in |36j . 

The parametrization in terms of Rf, the number of boxes in each column, naturally 
relates the crystal to the branes in section [3j This is the parametrization considered in 



29 Our considerations here are all in the closed string slicing as defined in [27] . There is also the open 
string slicing that is useful for relating (Wr) to the Wilson loop vev. 
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. It was found tha 





oo 




R M 




-Rm_i + 1 


(W R ) = (0| 


]I r +(^ 1 ) 








n 




.i=l 




i=i 




i=Rlj+2 



n r 



r + ( g -^+ M ); 



7V+M 



n r (9 



|0). (B.3* 



One can derive (I3.24p from this. 

The parametrization in terms of U is useful for relating the crystal to the bubbling 
geometry. (Wr) is written as 



(W R ) = (0| 

'hm-1 



X 



X 



n r +(^ 

i=l 

Y[ r_(g- ,i -' a --- J 



j=l 



2m-2 — 



. 1=1 

hm+1 



■'2 — •■■ f2m-l-'2m-! 



J Li=l 

'2m 

|0>. 



i=l 



i=l 



(B.39) 



Using ( IB. 361) . one can show that this reduces to (14.291) . The crystal is directly related to 
the bubbling geometry as one sees by comparing Figures 0] and [6j 



C Target space theory of open topological strings 

Let us consider the case of a single A-brane wrapping a Lagrangian submanifold L of 
a Calabi-Yau manifold M. In A-model, only (the cohomology class of) the symplectic 
structure u = ^u^dx 11 A dx v of M, not the whole Kahler structure, is relevant. When 
the B-field vanishes, the equation of motion for the gauge field is simply F = dA = 0. 
When the B-field is turned on the equation of motion should be generalized so that it 
is invariant under the gauge transformation B — » B + dA, A — > A — AFH. The minimal 
gauge invariant equation of motion is 



F + B = 0. 



(C.40) 



30 The expressions in (|B.37IIB.39|) differ in the parametrization of the Young diagram as well as the 
overall rescaling in the arguments of T±. Neither affects the amplitude. The overall rescaling is equivalent 
to conjugating all the operators by a power of q Lo , where Lq is the Virasoro zero mode |26) . 

31 The notation for gauge helds in this appendix differs from the rest of the paper. Here A is a hermitian 



connection, and A 



iA is its complexihcation. 
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For a single Lagrangian A-brane, one real modulus is given, when B = 0, by the 
holonomy § A along a non-trivial 1-cycle of L. This is again to be promoted to a gauge- 
invariant combination. The minimal combination is 



A+ B, (C.41) 

3D J D 

where D C M is a disk (or more generally a surface with boundary) whose boundary is 
the 1- cycled. 

In topological field theory, all moduli are complex. There must be another real modulus 
associated with the 1-cycle that combines with (1C.41I) to form a complex modulus. The 
existence of the additional modulus can be understood as follows [37] . 

Let us consider an infinitesimal deformation of L that preserves the condition that the 
submanifold is Lagrangian. An infinitesimal deformation corresponds to a normal vector 
field v, i.e., a section of the normal bundle TM\l/TL. We can think of v as a vector field 
defined in a neighborhood of L 

v = (C.42) 

(indices are for M coordinates), up to a vector field £, defined in the same domain, such 
that G TL. uj is preserved by the deformation if and only if 

L v to = (C.43) 

where L v is the Lie derivative. By contracting v with the symplectic form u, one can 
form a 1-form in the neighborhood of L: 

i v uj = uj flu v tM dx u . (C.44) 

i v uj is closed and its restriction i v u\ L defines a cohomology class of L. i v oj\l is exact if and 
only if v is a Hamiltonian vector field: v = {uj~ l )^ u d v f for some function /. Indeed it is 
known that two Lagrangian submanifolds related by a Hamiltonian deformation defines 
the same open string theory [37] . Infinitesimally, the deformation modulus is given by 

§ d D 

The modulus for a finite deformation is known to be given (for a single brane) by the 
symplectic area of the dislJ^I 

u. (C.45) 

D 



32 Since b\(M) = 0, the 1-cycle of L has to be the boundary of some chain D in M. 

33 When D is a holomorphically embedded disk, the symplectic area agrees with the area with respect 
to the Kahler metric of M. 
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One can easily check that the variation of (1Q.45j) is proportional to § dD i v uo. 

The description of the deformation modulus in the above two paragraphs is rather 
unsatisfactory because unlike (1C.41I) . ( ]C.45j) does not depend on a field of the world- 
volume theory. We now point out that the correct description of the brane, namely the 
Chern-Simons theory with a complex gauge group, generalizes (1C45|) into a form that is 
in parallel with ( 1C.41I) . 

In string theory, all deformations of the theory must appear in the physical spectrum. 
As Witten showed in [19], the spectrum comes solely from H l (L). For N coincident 
branes, this gives rise to U(N) theory when a reality condition is imposed on the string 
field. To account for deformations of the Lagrangian submanifold, however, it is more 
natural not to impose the reality condition [37] . Thus the gauge field is a complex field, 
and the gauge group is GL(N, C). We write A + i<p for the complex gauge field. 

To understand the meaning of the imaginary part of the gauge field, let us consider the 
vertex operators for the gauge fields. Let us use coordinates X a to parameterize L, and 
Y a for transverse directions so that X a + iY a are complex coordinates of M. Topological 
string theory is a version of bosonic string theory [19J. The "fixed" vertex operators [38] 
for the real gauge fields A a are the fermions x a m the notation of [19] . The corresponding 
"integrated" vertex operators are {bo, x a }, which turn out to be d T X a . (r is the parameter 
for the boundary of a world-sheet, and X a (z = a + it) are the embedding coordinates.) 
Thus the "integrated" vertex operators for the imaginary gauge fields cj) a are idtX a . In 
topological A-model, only the field configurations such that X a + iY a is holomorphic in 
z contribute. For such configurations id T X a = —id a Y a . It is well-known that integrating 
these vertex operators induces transverse deformations of the brane. Thus we have shown 
that the imaginary gauge fields give rise to deformations of the brane. 

Our conclusion, which is in agreement with some of the literature [37], is that in 
topological A-model the natural gauge group for the Chern-Simons theory on A-branes 
is a complex gauge group. This obviously raises the issue of non-perturbative definition 
because the path- integral seems to be divergent. Our attitude here is that only the 
perturbative definition of the Chern-Simons theory is relevant. This is consistent with 
the fact that topological string theory is defined only perturbatively. 

Once we accept the complex gauge group, there is natural generalization of (1C.45j) for 
the real modulus: 

I [ u (C.46) 

J 3D JD 

In other words, the complex modulus corresponding to the 1-cycle is 



A + J, (C.47) 

8D J D 

where A = <f) — iA, J = (u — iB). Note that this expression is invariant under the 
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complex generalization of the B-field gauge transformation 

B + iuj -> B + iuj + d(A + A + -+ A + — (A + iX). (C.48) 

The minimal non-Abelian gauge invariant equation of motion for the open strings 
interacting with closed strings is 

F + J1nxn = 0, (C.49) 

where T = dA + A 2 . 

For a promising approach to the target space theory, see [39J. 
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